BACKGROUND OF THE INVENTION 
1 . Field of the Invention 

The present invention relates to a process and a device 
for displacing a moveable unit on a base. 
5 Said device is of the type comprising a 

controllable actuator, for example an electric motor, 
intended to give rise to a linear displacement of the 
moveable unit on the base, as well as a system which is 
formed of a plurality of elements which are brought 
10 into motion upon the displacement of said moveable 
unit . 

Within the context of the present invention, 
said system exhibits at least two different motions and 
comprises as elements which may be brought into motion, 
15 in particular: 

- said base which can be mounted elastically with 
respect to the floor, especially so as to isolate it 
from vibrations originating from said floor; and/or 

- one or more auxiliary masses, for example measurement 
20 supports and/or loads, which are tied elastically to 

the base; and/ or 
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- one or more auxiliary masses, for example likewise 
measurement supports and/or loads, which are tied 
elastically to the moveable unit. 

When the moveable unit is set into motion, 
said elements of the system begin to move. However, 
especially by reason of the aforesaid elastic link, 
these elements still continue to move when the 
displacement of the moveable unit has terminated and 
when the latter comes to a stop. 

Such a continuance of the motions of said 
system is generally undesirable, since it may entail 
numerous drawbacks. In particular, it may disturb 
measurements, especially positioning measurements, 
which are made on the moveable unit or on these 
elements . 

Also, an object of the present invention is to 
control the moveable unit in such a way that all the 
moving elements of said system, for example the base 
and/or auxiliary masses, are stationary at the end of 
the displacement of the moveable unit. 

As regards said base, if it is mounted 
elastically with respect to the floor, it is known 

that, when the moveable unit is set into motion, during 
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the acceleration and deceleration phases, it is 
subjected to the reaction of the force applied to the 
moveable unit by the actuator. This reaction load 
excites the base which then oscillates on its supports. 
5 This disturbs the relative positioning of the moveable 
unit with respect to the base, and greatly impedes the 
accuracy of the device. 

This relative position error persists after 
the end of the displacement of the moveable unit and 
10 disappears only after the stabilization (which takes 
place much later) of the base. 

Various solutions for remedying this drawback 
are known. Some of these solutions make provision in 
particular : 

15 - to immobilize the base during the acceleration and 
deceleration phases via a disabling system, for 
example an electromagnetic disabling system, which is 
mounted in parallel with the elastic supports. 
However, this known solution prevents the supports 

20 from isolating the base from the vibrations 

originating from the floor during said ■ acceleration 
and deceleration phases ; 
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- to cancel the effect produced by the force developed 
by the actuator, by making provision for an 
additional actuator which is arranged between the 
base and the floor and which develops an additional 

5 force of the same amplitude but oppositely directed; 
or 

- to displace an additional moveable unit on the base 
according to a similar displacement, but oppositely 
directed, with respect to the displacement of the 

10 moveable unit, so as to cancel the inertia effects. 

However, none of these known solutions is 
satisfactory, since their effectivenesses are 
restricted and since they all require supplementary 
means (disabling system, additional actuator, 

15 additional moveable unit) which increase in particular 
the complexity, the cost and the bulkiness of the 
device . 

Moreover, above all, these solutions implement 
an action which acts only on the base and not on the 
20 other elements of the system which, for their part, 

continue to move when the moveable unit is stationary. 

The object of the present invention is to 

remedy these drawbacks. It relates to a process for 
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displacing, in an extremely accurate manner and at 
restricted cost, a moveable unit on a base mounted for 
example on the floor, whilst bringing all the motions 
to which this displacement gives rise to a stop at the 
end of the displacement, said moveable unit being 
displaced linearly according to a displacement which is 
predetermined in terms of distance and time, under the 
action of a controllable force. 

Accordingly, said process is noteworthy 
according to the invention in that: 

a) equations are defined which: 

- illustrate a dynamic model of a system formed by 
elements, of which said moveable unit is one, 
which are brought into motion upon a displacement 
of said moveable unit; and 

- comprise at least two variables, of which the 
position of said moveable unit is one; 

b) all the variables of this system, together with said 

force, are expressed as a function of one and the 

same intermediate variable y and of a specified 

number of derivatives as a function of time of this 

intermediate variable, said force being such that, 

applied to said moveable unit, it displaces the 
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latter according to said specified displacement and 
renders all the elements of said system immobile at 
the end of said displacement; 

c) the initial and final conditions of all said 
variables are determined; 

d) the value as a function of time of said intermediate 
variable is determined from the expressions for the 
variables defined in step b) and said initial and 
final conditions; 

e) the value as a function of time of said force is 
calculated from the expression for the force, 
defined in step b) and said value of the 
intermediate variable, determined in step d) ; and 

f) the value thus calculated of said force is applied 
to said moveable unit. 

Thus, the force applied to the moveable unit 

enables the latter to carry out the predetermined 

displacement envisaged, especially in terms of time and 

distance, whilst rendering the elements brought into 

motion by this displacement immobile at the end of the 

displacement so that they do not oscillate and, in 

particular, do not disturb the relative positioning 

between themselves and the moveable unit. 
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It will be noted moreover that, by reason of 
this combined control of said moveable unit and of said 
moving elements, one obtains an extremely accurate 
displacement of the moveable unit in a reference frame 
independent of the base and tied for example to the 
floor . 

It will be noted that the implementation of 
the process in accordance with the invention is not 
limited to a displacement along a single axis, but can 
also be applied to displacements along several axes 
which can be regarded as independent. 

Advantageously, in step a) , the following 
operations are carried out: the variables of the system 
are denoted xi, i going from 1 to p, p being an integer 
greater than or equal to 2, and the balance of the 
forces and of the moments is expressed, approximating 
to first order if necessary, in the so-called 
polynomial matrix form: 
A(s)X = bF 

with : 

• A(s) matrix of size p x p whose elements Aij (s) are 
polynomials of the variable s = d/dt; 



X the vector 



v X P , 
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• b the vector of dimension p; and 

• F the force exerted by the motor. 

Advantageously, in step b) , the following operations 
are carried out: 

- the different variables xi of said system, i going 
from 1 to p, each being required to satisfy a first 
expression of the form: 

j = r 

xi = I pi, j.y 01 , 

the y {3) being the derivatives of order j of the 
intermediate variable y, r being a predetermined 
integer and the pi, j being parameters to be 
determined, a second expression is obtained by 
putting y (:,) =s j .y: 

xi = ( E pi ' j - sj ) Y = pi(s> -y ' 

- a third expression of vectorial type is defined on 
the basis of the second expressions relating to the 
different variables xi of the system: 

X = P.y 



comprising the vector P = 



- said vector P is calculated, by replacing X by the 
value P.y in the following system: 
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B T .A(s) .P(s) = Op - 1 
< j=p 

bp . F = 2 Ap, j(s) .P j(s) .y 

j=i 

in which: 

. B T is the transpose of a matrix B of size px(p-l) , 
such that B T b=Op-l; 

. bp is the p-th component of the vector b previously 
defined; and 

. Op-1 is a zero vector of dimension (p-1) ; 

- the values of the different parameters pi,j are 
deduced from the value thus calculated of the vector 
P; and 

- from these latter values are deduced the values of 
the variables xi as a function of the intermediate 
variable y and of its derivatives, on each occasion 
using the corresponding first expression. 

Thus, a fast and general method of calculation 
is obtained for calculating the relations between the 
variables of the system and said intermediate variable, 
in the form of linear combinations of the latter and of 
its derivatives with respect to time. 

Advantageously, in step d) , a polynomial 
expression for the intermediate variable y is used to 
determine the value of the latter. 
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In this case, preferably, the initial and 
final conditions of the different variables of the 
system, together with the expressions defined in step 
b) , are used to determine the parameters of this 
5 polynomial expression. 

In a first embodiment, for displacing a 
moveable unit on a base which is mounted elastically 
with respect to the floor and which may be subjected to 
linear and angular motions, advantageously, the 
10 variables of the system are the linear position x of 
the moveable unit, the linear position xB of the base 
and the angular position 0z of the base, which satisfy 
the relations: 



x = y + 
xB = 
9z 
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15 in which: 

- m is the mass of the moveable unit; 

- mB, kB, k0, rB, r0 are respectively the mass, the 
linear stiffness, the torsional stiffness, the linear 
damping and the torsional damping of the base; 



12 

Atty Docket No. 155603-0209 ' 
Express Mail Label No. EL692571665US 



BJY/wrj 



- J is the inertia of the base with respect to a 
vertical axis; 

- d is the distance between the axis of translation of 
the center of mass of the moveable unit and that of the 
base; and 

- Y n) , Y {2) r y {3) and y {4) are respectively the first to 
fourth derivatives of the variable y. 



remedy 1 the aforesaid drawbacks (inaccurate 
displacement, etc) related to the setting of the base 
into oscillation during the displacement of the 
moveable unit. 



base a moveable unit on which are elastically mounted a 
number p of auxiliary masses MAi, p being greater than 
or equal to 1, i going from 1 to p, advantageously, the 
variables of the system are the position x of the moveable 
unit and the (linear) positions zi of the p auxiliary 
masses MAi, which satisfy the relations: 



This first embodiment makes it possible to 



In a second embodiment, for displacing on a 




— +— s + 1 
kj kj 
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in which: 

- IT illustrates the product of the associated 
expressions; 

- mi, zi, ki and ri are respectively the mass, the 
position, the stiffness and the damping of an auxiliary 
mass MAi; 

- m j , kj and rj are respectively the mass, the 
stiffness and the damping of an auxiliary mass MA j ; and 

- s=d/clt. 

In a third embodiment, for displacing a 
moveable unit on a base which is mounted elastically 
with respect to the floor and on which is elastically 
mounted an auxiliary mass, advantageously, the 
variables of the system are the positions x, xB and zA 
respectively of the moveable unit, of the base and of 
the auxiliary mass, which satisfy the relations: 

x = [(m As 2 + rAs + kA) . (mBs 2 + (rA + rB)s + (IcA + kB)) - (rAs + kA) 2 ].y 
< xB = -My (2) 

zA = -M(rAy (3 > +kAy (2) ) 
in which: 

- M, mB and mA are the masses respectively of the 
moveable unit, of the base and of the auxiliary mass; 

- rA and rB are the dampings respectively of the 

auxiliary mass and of the base; 
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- kA and kB are the stiffnesses respectively of the 
auxiliary mass and of the base; and 

- s=d/dt. 

In a fourth embodiment, for displacing on a 
base mounted elastically with respect to the floor, a 
moveable unit on which is elastically mounted an 
auxiliary mass, advantageously, the variables of the 
system are the positions x, xB and zC respectively of 
the moveable unit, of the base and of the auxiliary 
mass, which satisfy the relations: 

x = [(mCs 2 + rCs + kC).(mBs 2 +rBs + kB)].y 
<xB = [(mCs 2 + rCs + kC).(Ms 2 +rCs + kC)-(rCs + kC) 2 ].y 

zC = (rCs + kC).(mBs 2 +rBs + kB).y 

in which: 

- M, mB and mC are the masses respectively of the 
moveable unit, of the base and of the auxiliary mass; 

- rB and rC are the dampings respectively of the base 
and of the auxiliary mass; 

- kB and kC are the stiffnesses respectively of the 
base and of the auxiliary mass; and 

- s=d/dt. 

The present invention also relates to a device 
of the type comprising: 
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- a base mounted directly or indirectly on the floor; 

- a moveable unit which may be displaced linearly on 
said base; and 

- a controllable actuator able to apply a force to said 
moveable unit with a view to its displacement on said 
base . 

According to the invention, said device is 
noteworthy in that it furthermore comprises means, for 
example a calculator: 

- which implement steps a) to e) of the aforesaid 
process, so as to calculate a force which, applied to 
said moveable unit, makes it possible to obtain the 
combined effect or control indicated above; and 

- which determine a control command and transmit it to 
said actuator so that it applies the force thus 
calculated to said moveable unit, during a 
displacement . 

Thus, over and above the aforesaid advantages, 
the device in accordance with the invention does not 
require any additional mechanical means, thereby 
reducing its cost and its bulkiness and simplifying its 
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embodiment, with respect to the known and aforesaid 
devices , 

The figures of the appended drawing will 
elucidate the manner in which the invention may be 
5 embodied. In these figures, identical references 
designate similar elements. 

Figures 1 and 2 respectively illustrate two 
different embodiments of the device in accordance with 
the invention. 
10 Figures 3 to 7 represent graphs which 

illustrate the variations over time of variables of the 
system, for a first embodiment of the device in 
accordance with the invention. 

Figures 8 to 13 represent graphs which 
15 illustrate the variations over time of variables of the 
system, for a second embodiment of the device in 
accordance with the invention. 

The device 1 in accordance with the invention 
and represented diagrammatically in figures 1 and 2, 
2 0 according to two different embodiments, is intended for 
displacing a moveable unit 4, for example a moveable 
carriage, on a base 2, in particular a test bench. 
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This device 1 can for example be applied to 
fast XY tables used in microelectronics, to machine 
tools, to conveyors, to robots, etc. 

In a known manner, said device 1 comprises, in 
5 addition to the base 2 and to the moveable unit 4: 

- supports 3, of known type, arranged between the base 
2 and the floor S; 

- means (not represented), for example a rail, fixed on 
the base 2 and enabling the moveable unit 4 to be 

10 displaced linearly on said base 2; and 

- a controllable actuator 5, preferably an electric 
motor, able to apply a force F to said moveable unit 
4 with a view to its displacement on the base 2. 

Within the context of the present invention, 
15 the device 1 comprises a system SI, S2 which is formed 
of various elements specified hereinbelow and variables 
according to the embodiment contemplated, which are 
brought into motion upon the displacement of the 
moveable unit 4 . 
20 According to the invention, said device 1 is 

improved in such a way as to obtain directly at the end 
of a displacement of the moveable unit 4: 
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- accurate positioning of the latter in a reference 
frame (not represented) , independent of the moveable 
unit 4 and of the base 2 and tied for example to the 
floor; and 

5 - immobilization of all the moving elements of said 
system SI , S2 . 

To do this, the device 1 moreover comprises, 
according to the invention, calculation means 6 which 
calculate a particular force F, which is intended to be 
10 transmitted in the form of a control command to the 
actuator 5, as illustrated by a link 7, and which is 
such that, applied to said moveable unit 4, it produces 
a combined effect (and hence combined control) : 

- on the one hand, on the moveable unit 4 so that it 
15 exactly carries out the envisaged displacement, 

especially as regards the prescribed duration and 
prescribed distance of displacement; and 

- on the other hand, on said system Si, S2 so that all 
its moving elements are immobile at the end of the 

20 displacement of the moveable unit 4. 
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Accordingly, said calculation means 6 
implement the process in accordance with the invention, 
according to which: 

a) equations are defined which: 

5 - illustrate a dynamic model of said system (for 

example Si or S2) formed by the different 
elements, of which said moveable unit 4 is one, 
which are brought into motion upon a displacement 
of said moveable unit 4; and 
10 - comprise at least three variables, of which the 

position of said moveable unit 4 is one; 

b) all the variables of this system, together with said 
force F, are expressed as a function of one and the 
same intermediate variable y and of a specified 

15 number of derivatives as a function of time of this 
intermediate variable, said force F being required 
to be such that, applied to said moveable unit 4, it 
displaces the latter according to said specified 
displacement and renders all the elements of said 

2 0 system immobile at the end of said displacement; 

c) the initial and final conditions of all said 
variables are determined; 
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d) the value as a function of time of said intermediate 
variable is determined from the expressions for the 
variables defined in step b) and said initial and 
final conditions; and 
5 e) the value of said force is calculated from the 

expression for the force, defined in step b) and 
said value of the intermediate variable, determined 
in step d) . 

Thus, by virtue of the invention, the force F 
10 applied to the moveable unit 4 enables the latter to 
carry out the predetermined displacement envisaged, 
especially in terms of time and distance, whilst 
rendering the elements (specified hereinbelow) which 
are brought into motion by this displacement immobile 
15 at the end of the displacement so that they do not 
oscillate and, in particular, do not disturb the 
relative positioning between themselves and the 
moveable unit 4. 

It will be noted moreover that, by reason of 
20 this combined effect or control of said moveable unit 4 
and of said moving elements, one obtains an extremely 
accurate displacement of the moveable unit 4 in a 
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reference frame independent of the base 2 and tied for 

example to the floor S. 

Of course, the implementation of the present 

invention is not limited to a displacement along a 
5 single axis, but can also be applied to displacements 

along several axes which can be regarded as 

independent . 

According to the invention, in step d) , a 

polynomial expression for the intermediate variable y 
10 is used to determine the value of the latter, and the 

initial and final conditions of the different variables 

of the system, together with the expressions defined in 

step b) are used to determine the parameters of this 

polynomial expression . 
15 The process in accordance with the invention 

will now be described in respect of four different 

systems (of moving elements) . 

In a first embodiment (not represented) , the 

supports 3 are of elastic type and make it possible to 
20 isolate the base 2 from the vibrations originating from 

said floor S. The natural frequency of the base 2 on 

said elastic supports 3 is generally a few Hertz. 

Furthermore, in addition to the translational motion of 
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the moveable unit 4 controlled by the force F, an 
angular motion is created between the base 2 and the 
moveable unit 4. Specifically, in this case, the axis 
of the moveable unit 4 does not pass through its center 
5 of mass, the force produced by the actuator 5 creates a 
moment about the vertical axis. The rail is assumed to 
be slightly flexible and thus allows the moveable unit 
4 small rotational motions about the vertical axis, 
which corresponds to the aforesaid relative angular 
10 motion between the base 2 and the moveable unit 4. 



displace the moveable unit 4 on the base 2 which is 
mounted elastically with respect to the floor and which 
may be subjected to a (relative) angular motion, the 
15 variables of the system are the linear position x of 

the moveable unit 4, the linear position xB of the base 
2 and the angular position 0z of the base 2, which 
satisfy the relations: 



Consequently, in this first embodiment, to 




kBke kBke 



rBJ mBre 




20 



in which 
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- m is the mass of the moveable unit 4; 

- mB, kB, k9, rB, r9 are respectively the mass, the 
linear stiffness, the torsional stiffness, the linear 
damping and the torsional damping of the base 2; 

5 - J is the inertia of the base 2 with respect to a 
vertical axis; 

- d is the distance between the axis of translation of 
the center of mass of the moveable unit 4 and that of 
the base 2; and 

10 - y (1) , y {2) , y (3) and y (4) are respectively the first to 
fourth derivatives of the variable y. 

Specifically, in this first embodiment, the 

balance of the forces and of the moments, the angle Qz 
being approximated to first order, may be written: 

mx< 2) = F 

15 \ mBxB (2) =-F-kBxB-rBxB (1) (1) 

J9z (2) =-dF-keez-reez (1) 

It will be noted that, within the context of 

the present invention, a (p) is the derivative of order (3 
with respect to time of the parameter a, regardless of 
a. Thus, for example, x (1) is the first derivative of x 
20 with respect to time. 
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The calculation of the intermediate variable y 



is achieved by putting s = — , x = P(s)y, xB = PB(s)y, 

dt 

Oz = PG (s)y and by rewriting the system (1) with this 
notation: 

ms 2 P(s)y = F 

(mBs 2 +rBs + kB)PB(s)y = -F 
(Js 2 +r0s + k9)P9(s)y = -dF 



i.e.: 



(mBs 2 +rBs + kB)PB(s) = - (Js 2 + rGs + k9) P6 (s) = -ms 2 P(s) 



and hence: 
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From these expressions, we immediately deduce 
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The expression for y as a function of x, x (1) , 

xB, xB ,1) , 0z and 6z {1) is obtained by inversion. However, 

15 this formula is not necessary in order to plan the 
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trajectories of x, xB and Qz . Specifically, since we 
want a stop-stop displacement of the moveable unit 4 
between xO at the instant tO and xl at the instant tl, 
with 

5 x (1) (tO) = 0 = x (1) (tl) and xB(tO) = 0 = xB(tl) , xB (1) (tO) = 0 = xB (1) {tl) an 
d 0z(tO) = 0 = 0z(tl) , 9z (1) (t0) = 0 = 9z (1) (tl) , with in addition 
F(tO)=0=F(tl) , 

we deduce therefrom through the aforesaid expressions 
(2) that y(tO)=xO,y(tl)=xl and 
10 y a) (ti)=y (2, (ti) - y (3) (ti) - y (4) (ti) - y (5) (ti) = y (6) (ti) = 0, i = 0.1 
i.e. 14 initial and final conditions. 

It is sufficient to choose y as a polynomial 

with respect to time of the form: 

P 

y(t) - xO + (xl - xO) (G(t) ) a 2 a i (a(t) ) x (3 ) 

i = 0 

15 with 

t - to 



a(t) = 



tl - to 

and a > 7 and |3 > 6. The coefficients aO, . . . , a(3 are 
then obtained, according to standard methods, by- 
solving a linear system. 



26 

Atty Docket No. 155603-0209 

Express Mail Label No. EL692571665US 



BJY/wrj 



The reference trajectory sought for the 
displacement of the moveable unit 4 is then given by 
expressions (2) with y(t) given by expression (3). 

Moreover, the force F as a function of time to 
5 be applied to the means 5 is obtained by integrating 
the value of y obtained via expression (3) in the 

expression F(t) = M.x (2) (t) . 

In this first embodiment, we obtain: 



F(t) = M 



Y 2 + 



rB_ rQ_ 

kB + k9 



\ 



y< 3 > + 



f mB rBrG J ^ 



• + + ■ 



v kB kBkG k9 , 



y (4) + 



rBJ + mBr0^ 



v kBkG j 



(5) mBJ {6) 
kBkG 



(3A) 



10 with y(t) given by expression (3) . 

Thus, since by virtue of the device 1 the base 
2 is immobilized at the end of the displacement, it 
does not disturb the positioning of the moveable unit 4 
in the aforesaid reference frame so that said moveable 

15 unit 4 is positioned in a stable manner as soon as its 
displacement ends. Moreover, since its displacement is 
carried out in an accurate manner, its positioning 
corresponds exactly in said reference frame to the 
sought-after positioning. 

2 0 Represented in figures 3 to 7 are the values 

respectively of said variables y (in meters m) , x (in 
meters m) , xB (in meters m) , 0z (in radians rd) and F 
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(in Newtons N) as a function of time t (in seconds s) 
for a particular exemplary embodiment, for which: 

- m=40 kg; 

- mB=800 kg; 

5 - kB=mB (5 . 2n) 2 corresponding to a natural frequency of 
5 Hz; 



- rB=0 . 3VkBmB corresponding to a normalized damping of 
0.3; 

- J=12 0 Nm corresponding to the inertia of the moveable 
10 unit 4; 

- kQ=J (10 . 271) 2 corresponding to a natural rotational 
frequency of 10 Hz; 

- r »=0 . 3 Vk0J corresponding to a normalized rotational 
damping of 0.3; 

fl 15 - d=0.01 m corresponding to the of f -centering of the 
moveable unit 4; 

- tl-tO = 0.4 s; and 

- xl-xO = 25 mm. 
The moveable unit 4 is displaced from the 

20 position xO at rest (x0 (1) -0) at the instant t0, to the 

position xl at rest (xl (1) =0) at the instant tl . It is 

therefore displaced over a distance of 25 mm in 0.4 s. 
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To obtain this displacement, as well as the 
immobilization (at the end of said displacement) of the 
various motions to which the displacement gives rise, 
the force F represented in figure 7 should be applied 
5 to said moveable unit 4. This force is given by 

expression (3A) with y given by (3) for a = 7 and 
P = 6. In this case, the coefficients aO up to a6 are 
given by a0=1716, al=-9009, a2=20020, a3=-24024, 
a4=16380, a5=-6006, a6=924. 

10 In a second embodiment represented in figure 

1, the system SI comprises, in addition to the moveable 
unit 4, a number p of auxiliary masses MAi, p being 
greater than or equal to 1, i going from 1 to p, which 
are linked respectively by elastic links el to ep of 

15 standard type, in particular springs, to said moveable 
unit 4. In the example represented, p=3 . 

In this case, the variables of the system are 
the position x of the moveable unit 4 and the positions 
zi of the p auxiliary masses MAi, which satisfy the 

20 relations: 



Atty Docket No. 155603-0209 

Express Mail Label No. EL692571665US 



29 



BJY/wrj 



X = 



f p 'mi 2 ri ^ 
— s z + — s + 1 



n 

i=1 



ki 



ki 



zi = 



p y mj 2 rj ^ 

— +— s + 1 



n 

j=1 



kj kj 



■y 



— s + 1 
ki 



(4) 



in which: 



- IT illustrates the product of the associated 



expressions ; 

5 - mi, zi, ki and ri are respectively the mass, the 

position, the stiffness and the damping of an auxiliary 
mass MAi; 

- m j , kj and rj are respectively the mass, the 
stiffness and the damping of an auxiliary mass MA j ; and 



10 - s=d/dt. 



Specifically, the dynamic model of the system 



SI may be written: 



Mx (2) =F + £(ki(zi-x) + ri(zi (1) -x {1) )) 

i=1 

Mizi {2) = ki (x - zi) + ri (x (1) - zi {1) ), i = 1, . . . p. 



(5) 



As in the foregoing, we wish to find laws of 
15 motion which ensure the desired displacement of the 

moveable unit 4, the auxiliary masses MAi (for example 
measurement devices and/ or loads) being immobilized as 



soon as the moveable unit 4 stops 
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10 



Accordingly, the intermediate variable y is 
calculated by the same approach as earlier and the 
trajectory of the moveable unit 4 is planned by way 
thereof . 

The intermediate variable y being required to 



satisfy x=P(s)y, zi=Pi(s)y, i=l, . 



. p, with s = — , we 

dt 



must have, substituting these relations into the system 
(5) : 

(mis 2 +ris+ki) Pi = (ris+ki)P, i = l, . p 

From this expression, we immediately derive: 



P(s) = 



n 

i=1 



mi 2 n 

— s + — S + 1 

ki ki 











, Pi = 




n 







mi 2 n A 

— s z + — s + 1 
kj kj 



— S + 1 

ki 



thereby proving the aforesaid formulae (4) . 

In this case, it may be demonstrated that the 
force F to be applied satisfies the relation: 



15 F(t) = 



r r 
Ms 2 +' 



— s z + — s + 1 



kj kj 



The aforesaid formulae are verified and 
specified hereinbelow for two and three auxiliary 
masses MAi respectively. 
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In the case of two auxiliary masses (p=2) , the 



model may be written: 

Mx (2) = F - k1(x - z1) - r1(x {1) - z1 (1) ) - k2(x - z2) - r2(x (1) - z2 (1) ) 
m1z1 (2) = k1 (x - z1) + r1 (x (1) - z1 (1) ) 
m2z2 (2) = k2 (x - z2) + r2 (x (1) - z2 (1) ) 

From this we immediately deduce: 



x = ( — s z + — s + 1) ( — s z + — s + 1)y 



k1 



k1 



k2 



k2 



r1 .m2 2 r2 



z1 = (— s + 1) (— s z + — s + 1)y 
k1 k2 k2 

_ ,r2 „. ,m1 2 r1 
z2= — s + 1 — s + — s + 1y 
k2 k1 k1 



(6) 



nil ,o ri „ . . 
i.e, putting — = Ti^ and — = 2DiTi, 

ki ki 



i = 1.2 : 



x = y + 2(D1T1 + D2T2)y (1) + (Ti 2 + T2 2 + 4D1D2T1T2)y (2) 
+ 2(D1T1T2 2 + D2T2T1 2 )y (3) + (T1 2 T2 2 )y l4) 

1 = y + 2(D1T1 + D2T2)y (1) +(T2 2 + 4D1D2T1T2)y (2) + (2D1T1T2 2 )y <3) 

2 = y + 2(D1T1 + D2T2)y (1) + (T1 2 +4D1D2T1T2)y (2) +(2D2T2T1 2 )y (3) . 



The expression for y, or more precisely the 



expressions for y, y (1) , y <2) , y (3) , y (4) and y V3 ' , are 
10 deduced therefrom by inverting the system obtained on 



.(5) 



the basis of x, zl, z2, x (1) , zl li, ( z2 l± \ 

We deduce therefrom that, to perform a 
displacement from xO at the instant tO to xl at the 
instant tl, with the auxiliary masses at rest at tO and 
15 tl, it is sufficient to construct a reference 



(i) 



(1! 
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trajectory for y with the initial and final conditions 
y(tO)=xO, y(tl)=xl and all the derivatives 
y (k) (tO) =y (k) (tl) =0, k varying from 1 to 6 or more if 
necessary, and to deduce therefrom the reference 
trajectories of the main and auxiliary masses, as well 
as of the force F to be applied to the motor. 

In this case, the force F satisfies the 

relation: 



Furthermore, the model for three auxiliary 
masses MAi (p=3) [see figure 1], may be written, as 
earlier : 



Mx ,2) = F - k1(x - z1) - r1(x (1> - z1 (1) ) 

- k2(x - z2) - r2(x (1) - z2 (1) ) - k3(x - z3) - r3(x (1) - z3 (1) ) 
< m1z1 (2> = k1(x - z1) + r1(x (1) - z1 (1) ) 
m2z2 (2) = k2(x - z2) + r2(x n) - z2 (1) ) 
m3z3 (2) = k3(x - z3) + r3(x (1) - z3 (1) ) 




k2 k2 



— s + — s + 1 
{ k2 k2 J 

ml 2 r1 \ 
— s 2 +— s + 1 ]y. 

k1 kl 



From this we immediately deduce: 
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x = ( — + — s + 1)( — s z + — s + 1)(— s* + — s + 1)y 



k1 

z1 = (lls + 1)( 
k1 



k1 k2 

m2 2 ""2 ... 
— + — s + 1 
k2 k2 



,r2 „,,m1 2 r1 



k2 k3 

m3 2 ""3 
— s + — 
k3 k3 

m3 2 r3 



k3 
s + 1)y 



(7) 



z2 = (— s + 1)( — s z + — s + 1)( s z + — s + 1)y 

k2 k1 k1 k3 k3 

_ ,r3 „ w m1 2 r1 ..w m2 2 , f2 , 14 
z3= — s + 1 — s^ + — s + 1 — s z + — -s + 1)y 
k3 k1 k1 k2 k2 



We proceed as earlier in order to determine 
the values as a function of time of the different 
variables and in particular of the force F, the latter 



5 satisfying the expression: 

F(t) = [(Ms 2 + (r1 + r2 + r3)s + (k1 + k2 + k3)). 



ml 
k1 



s + — s + 1 



k1 



s 2 + — s + 1 



- (r1s + kl) 



m2 



k2 
r2 



k2 



m3 2 r3 

s z + — s + 1 



k3 



•s' + — s + 1 
v k2 k2 

(r2s + k2)f^s 2 +-s + 1 
v \ k1 k1 

^m! 2 r1 
— s +— s + 1 
k1 k1 



Ym3 2 r3 



k3 



j 



s + — s + 1 
A k3 k3 

^^m3 2 r3 



- (r3s + k3) 



< 



s^ + — s + 1 
k3 k3 

m2 2 r2 
— s + — s + 1 
k2 k2 



]y. 



Represented in figures 8 to 13 are the values 
reS p ec tively of the variables y, x, zl, z2 , z3 and F as 
a function of time t for a particular example of the 
10 embodiment of figure 1, zl to z3 being the 

displacements of the auxiliary masses MAl , MA2 and MA3 
respectively. The variables y, x, zl, z2 and z3 are 
expressed in meters (m) and the force F in Newtons (N) . 
This example is such that: 



15 - M=5 kg; 
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- ml = 0 . 1 kg; 

- m2=0.01 kg; 

- m3=0 . 5 kg; 

- kl=ml (5.271) 2 , k2=m2 (4.2ti) 2 , k3=m3 ( 6 . 2n) 2 , 

5 corresponding to natural frequencies of 5, 4 and 6 Hz 
respectively; 

- rl = 0.3Vklml , r2 = 0 . 2Vk2m2 , r3 = 0 . 15Vk3m3 , corresponding 
to normalized dampings of 0.3, 0.2 and 0.15 
respectively; 

10 - tl-tO = 0.34 s; and 

- xl-xO = 40 mm. 

Additionally, in a third embodiment 
represented in figure 2, the system S2 comprises the 
moveable unit 4, the base 2 which is mounted 

15 elastically with respect to the floor S and an 

auxiliary mass MA which is linked by way of an elastic 
link eA of standard type to said base 2. 

In this case, the variables of the system are 
the positions x, xB and zA of the moveable unit 4, of 

2 0 the base B and of the auxiliary mass MA, which satisfy 
the relations: 
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x- [(mAs 2 +rAs + kA).(mBs 2 + (rA + rB)s + (kA + kB)) - (rAs + kA) 2 ].y 
< xB = -My (2) (8) 
z A = -M (r Ay i3} +kAy (2) ) 

in which: 

- M, mB and mA are the masses respectively of the 
moveable unit 4, of the base 2 and of the auxiliary 

5 mass MA; 

- rA and rB are the dampings respectively of the 
auxiliary mass MA and of the base 2; 

q - kA and kB are the stiffnesses respectively of the 

M auxiliary mass MA and of the base 2; and 

yj 10 - s=d/dt . 

f;j Specifically, the dynamic model of the system 

H S2 may be written: 

E \ Mx (2) = F 

tl mBxB ,2) =-F-kBxB-rBxB (1) -k(xB-zA)-rA(xB (1) -zA ,1) ) (9) 

mzA (2) = kA(xB - zA) + rA(xB™ - zA (1) ) 

The intermediate variable must satisfy: 

d 

15 x = P(s)y, xB = PB(s)y and zA = Pz(s)y with s = — . 

dt 

Substituting these expressions into (9) , we 

obtain: 
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5 



F = Ms 2 P(s)y 

(mBs 2 + (rA + rB)s + (kA + kB)) PB(s) = -Ms 2 P(s) + (rAs + kA) Pz(s) 
(mAs 2 + rAs + kA) Pz(s) = (rAs + kA) PB(s). 

On eliminating Pz from the last equation, it 
follows that: 

[(mAs 2 + rAs + kA)(mBs 2 + (rA + rB)s + (£A + kB)) - (rAs + kA) 2 ]PB = ~(mAs 2 + rAs + kA)Ms 2 P 

from which we derive: 

P = (mAs 2 + rAs + kA)(mBs 2 + (rA + rB)s + (kA + kB)) - (rAs + kA) 2 
PB = -Ms 2 
Pz = -Ms 2 (rAs + kA) 

thus making it possible to obtain the aforesaid 
expressions (8) . 
10 The values as a function of time of the 

different variables, and in particular the force F, are 
then obtained as before. 

In this case, said force F satisfies the 
expression: 

1 5 F(t) - M [ (mAs 2 + rAs + kA)(mBs 2 + (rA + rB)s + (kA + kB))- (rAs + kA) 2 ] y <2) . 

In a fourth and last embodiment (not 
represented) , the system is formed of the moveable unit 
4, of the base 2 and of an auxiliary mass MC which is 
tied elastically to said moveable unit 4. 
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In this case, the variables of the system are 



the positions x, xB and zC respectively of the moveable 
unit 4, of the base 2 and of the auxiliary mass MC, 
which satisfy the relations: 



zC = (rCs + kC)(mBs 2 + rBs + kB).y 

in which: 

- M, mB and mC are the masses respectively of the 
moveable unit 4, of the base 2 and of the auxiliary 
mass MC; 

10 - rB and rC are the dampings respectively of the base 2 
and of the auxiliary mass MC; 

- kB and kC are the stiffnesses respectively of the 
base 2 and of the auxiliary mass MC; and 

- s=d/dt. 

15 Specifically, the dynamic model of this system 

may be written: 



5 



x = [(rnCs 2 + rCs + kCUmBs 2 + rBs + kB].y 
< xB = [ (rnCs 2 + rCs + kC) . (Ms 2 + rCs + kC) - (rCs + kC) 2 ].y 



Mx (2) = F - kC(x - zC) - rC(x (1 > - zC (1) ) 



mBxB U) = -F-kBxB-rBxB 



(10) 



mCzC (2) =kC(x-zC) + rC(x (1) -zC (1) ) 
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By using, as in the foregoing, the polynomial 
representation of the variable s ~"^' t -' ie s Y stem (10) 

becomes : 

(Ms 2 + rCs + kC)x = F + (rCs + kC)zC 
<(mBs 2 +rBs + kB)xB = -F 
(mCs 2 + rCs + kC)zC = (rCs + kC)x 

which, together with the expressions for each of the 
variables as a function of the intermediate variable 
(and of its derivatives), x = P(s)y, xB = PB(s)y, zC = Pz(s)y, 
finally gives: 

P - (mCs 2 + rCs + kC)(mBs 2 + rBs + kB) 
< PB = (mCs 2 + rCs + kC)(Ms 2 + rCs + kC) - (rCs + kC) 2 
Pz = (rCs + kC) (mBs 2 + rBs + kB) 

The construction of the reference trajectories 
of Yi and then of x, xB, zC and F is done as indicated 
earlier . 

In this case, the force F satisfies: 
F(t) = - (mBs 2 + rBs + kB)[ (mCs 2 + rCs + kc)(lVls 2 + rCs + kc)- (rCs + kC) 2 ] y. 

A method in accordance with the invention will 
now be described which makes it possible to determine 
in a general and fast manner the expressions defined in 
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the aforesaid step b) of the process in accordance with 
the invention, for linear systems of the form: 

£Ai,j(s)xj = biF, i = 1, ...,p (11) 

where the Ai,j(s) are polynomials of the 
5 variable s, which, in the case of coupled mechanical 
systems, are of degree less than or equal to 2 and 
where one at least of the coefficients bi is non-zero. 
F is the control input which, in the above examples, is 
the force produced by the actuator 5 . 
10 Accordingly, according to the invention, in 

step b) , the following operations are carried out: 

- the different variables xi of said system (for 
example Si or S2 ) , i going from 1 to p, p being an 
integer greater than or equal to 2 , each being 
15 required to satisfy a first expression of the form: 

xi = 2 pi, j.y , 

j=o 

the y <j) being the derivatives of order j of the 
intermediate variable y, r being a predetermined integer 
and the pi, j being parameters to be determined, a second 
20 expression is obtained by putting y (:j) =s j .y: 



xi = 



2 pi, j.s 3 

j=o 



y = Pi(s) .y 
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- a third expression of vectorial type is defined on the 
basis of the second expressions relating to the different 
variables xi of the system: 

X=P.y 



comprising the vectors 





f pl l 


p = 


>J 


X - 


. xp ; 



- said vector P is calculated, replacing X by the value P.y 
in the following expressions: 



B\A(s) .P(s) = Op - 1 

(3=P \ 

bp.F - S Ap, j(s) .P j(s) .y 

j=i 



in which: 

10 . B T is the transpose of a matrix B of size px(p-l) and of 
rank p-1, such that B T b=Op-l; 
* bp is the p-th component of the vector b; and 
. Op-1 is a zero vector of dimension (p-1) ; 

- the values of the different parameters pi,j are deduced 
15 from the value thus calculated of the vector P; and 

- from these latter values are deduced the values of the 
variables xi as a function of the intermediate variable y 
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and of its derivatives, on each occasion using the 



corresponding first expression. 

The aforesaid method is now justified. 

Let us denote by A(s) the matrix of size pxp whose 

5 coefficients are the polynomials Ai,j(s), i,j=l, p, 



i.e.: 



20 









f bl l 


, x - 




and b = 




J 









'a1,1(s) - Al,p(s) 
A(s) = \ i 

^Ap,l(s) Ap, p(s) 

Without loss of generality, it can be assumed that 
10 the rank of A(s) is equal to p (otherwise, the system is 
written together with its redundant equations and it is 
sufficient to eliminate the dependent equations) and that 
bp#0 . There then exists a matrix B of size px(p-l) and of 
rank p-1 such that: 
15 B T b=0p-l 

where T represents transposition and Op-1 the vector of 
dimension p-1, all of whose components are zero. The 
system (11) premultiplied by B T then becomes: 

B T A(s)X - Op - 1, bpF = £ Ap, jxj. 



(12) 



j=i 



As indicated earlier, an intermediate variable 

y is characterized in that all the components of the 

vector X can be expressed as a function of y and of a 
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finite number of its derivatives. For a controllable 
linear system, such an output always exists and the 
components of X can be found in the form of linear 
combinations of y and of its derivatives, i.e.: 



5 xi=£pi,jy 

i=o 



where y (j) is the derivative of order j of y with respect 
to time and where the pi,j are real numbers which are 
not all zero, or alternatively: 



10 



xi = ( j pi, js J )y = Pi(s)y, i = 1, . . . , p. 

]=0 



We shall calculate the vector P(s)= 



f P1(s) x 



v Pp(s) y 



, by 



replacing X by its value P(s)y in (12) : 



B T A(s)P(s) = Op - 1, bpF = J Ap, j(s)Pj(s)y. 



(13) 



j=i 



Consequently, P belongs to the kernel of the 



matrix B T A(s) of dimension 1, since B is of rank p-1 



15 and A(s) of rank p. To calculate P, let us denote by 
Al(s), Ap(s) the columns of the matrix A(s) and 

A(s) the matrix of size (p-l)x(p-l) defined by: 



A(s)=(A2(s) , . . . , Ap(s) ) . 
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Let us also denote by P(s) the vector of 



dimension p-1 defined by: 



P(s) = 



(P2[s) 



,Pp(s) y 



Let us rewrite (13) in the form 



5 B T A1 (s) PI (s) +B T A(s) P (s) =0p-l or alternatively 

B T A(s)P(s) = - B T Al(s)Pl(s) . Since the matrix B T A(s) is 
invertible, we have: 

P(s) = -(B T A(s)r 1 B T A1(s)P1(s) 



10 



i.e.: 
P(s) = - 



1 



det(B T A(s)) 



(co(B T A(s))) T B T A1(s)P1(s) 



(14) 



where co(B A(s)) is the matrix of the cof actors of B A ( s ) 

From this we immediately deduce that it is 
sufficient to choose: 



P1(s) = det(B T A(s)) 
P(s) = -(co(B T A(s))) T B T A1(s) 



(15) 



15 this completing the calculation of the vector P(s). 

It will be observed that if the Ai,j(s) are 
polynomials of degree less than or equal to m, the 
degree of each of the components of P is less than or 
equal to mp. Specifically, in this case, the degree of 
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the determinant det ( B T A ( s ) ) is less than or equal to 
(p-l)m and the degree of each of the rows of 
(co (B T A(s) ) ) T B T Al(s) 7 using the fact that the degree 
of a product of polynomials is less than or equal to 
5 the sum of the degrees, is less than or equal to (p- 
l)m+m=pm, hence the aforesaid result. 

In all the examples presented earlier, which 
model mechanical subsystems, we have m=2 . 

It may easily be verified that this general 
10 method yields the same calculations for P as in each of 
the examples already presented hereinabove. 

We shall return to certain of the examples 
dealt with earlier and show how the calculation of the 
variable y makes it possible to achieve passive 
15 isolation of the elastic modes. 

In all these examples, the trajectories are 
generated on the basis of polynomial trajectories of 
the intermediate value y, which are obtained through 
interpolation of the initial and final conditions. 
20 Furthermore, we are interested only in the particular 
case where the system is at rest at the initial and 
final instants, thereby making it possible to establish 
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simple and standard formulae which depend only on the 
degree of the polynomial. 

In the simplest case, where the initial and 
final derivatives of y are zero up to order 4, the 
5 sought-after polynomial is of degree 9: 

y(tO) - yO y {1) (tO) = 0 y {2) (tO) = 0 y (3) (tO) = 0 y (4) (tO) - 0 
y(t1) - y1 y (1, (t1) = 0 y (2} (t1) = 0 y< 3) (t1) = 0 y (4> (tl) = 0 

which gives : 



y(t) = yO + {y1 - y0)a 5 (1 26 - 420a + 540a 2 - 31 5a 3 + 70a 4 ), a = 



t-tO 
t1-t0 



(16) 



If we ask for a polynomial such that the 
10 initial and final derivatives are zero up to order 5, 
the sought-after polynomial is of degree 11: 

y(t) = yO + (y1-y0)a 6 (462-1 980a + 3465a 2 -3080a 3 + 1386a 4 -252a 5 ) 

still with a defined as in (16) . 

If we ask for a polynomial such that the 

15 initial and final derivatives are zero up to order 6, 

the sought-after polynomial is of degree 13: 

y (t) = yO + (y1 - y0)a 7 (1 7 1 6 - 9009a + 20020a 2 - 24024a 3 + 1 6380a 4 - 6006a 5 + 924a 6 ) . 
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